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Is it possible to increase T c by constructing cuprate heterostructures, which combine the high 
pairing energy of underdoped layers with the large carrier density of proximate overdoped layers? 
We investigate this question within a model bilayer system using an effective theory of the doped 
Mott insulator. Interestingly, the question hinges on the fundamental nature of the superconducting 
state in the underdoped regime. Within a plain slave boson mean field theory, there is absolutely no 
enhancement of T c . However, we do get a substantial enhancement for moderate inter-layer tunneling 
when we use an effective low energy theory of the bilayer in which the effective quasiparticle charge 
in the underdoped regime is taken as an independent phenomenological parameter. We study the 
T c enhancement as a function of the doping level and the inter-layer tunneling, and discuss possible 
connections to recent experiments by Yuli et al. [Phys. Rev. Lett. 101, 057005 (2008)]. Finally, we 
predict a unique paramagnetic reduction of the zero temperature phase stiffness of coupled layers, 
which depends on the difference in the current carried by quasiparticles on the two types of layers 
as (Ji — J2) 2 . 



PACS numbers: 74.78.Fk, 74.72. Dn, 74.20.Fg 

I. INTRODUCTION 

There are strong indications^^ that the supercon- 
ducting transition temperature of underdoped cuprate 
materials is limited only by their small supcrfluid density, 
while the pairing scale is very high. This understanding 
has raised the hope that T c can yet be made substan- 
tially higher by clever design of materials. In particular, 
a number of recent theoretical^^ and experimental^^ 
studies explored the possible benefit in heterostructured 
materials combining metallic layers with layers of under- 
doped cuprate material. The basic idea is simple; un- 
derdoped layers contribute a strong microscopic pairing 
interaction, whereas metallic layers provide a high den- 
sity of charge carriers. But it could also go the other 
way. Namely, the metal destroys pairing in the under- 
doped layer without contributing much of its charge car- 
riers. The question how much, if at all, such systems 
can actually enhance T c may require deeper knowledge of 
the nature of the superconducting state in the cuprates. 
Of particular importance in this respect is better under- 
standing of the mechanisms that reduce the superfluid 
density with temperature and with proximity to the Mott 
insulating state. 

In this paper we investigate the problem of T c enhance- 
ment within a bilayer model using an effective descrip- 
tion of the doped Mott insulator. The model captures 
an interesting competition of effects, which we expect is 
rather general to cuprate heterostructures and possibly 
inhomogeneous realizations of these materials. In par- 
ticular we discuss possible implications of our results to 
recent experiments in La 2 _xSr x Cu04 (LSCO) bilayer*^. 
We argue that experiments with heterostructures may 
shed new light on fundamental questions concerning the 
nature of superconductivity in cuprates. 

The essential physics that determines T c of the bilayer 




FIG. 1: Illustration of the two dimensional low energy dis- 
persion of a bilayer consisting of a d-wave superconducting 
layer and a nominally metallic layer. The four red surfaces 
mark the original Dirac cones of the superconducting layer. 
The blue surfaces are the Dirac cones induced by the prox- 
imity effect on the metallic layer. The blue curve denotes the 
original Fermi surface of that layer. 



is most clearly illustrated within the effective low energy 
theory of d-wave superconductor s] 11 ! 12 ' 13 . The superfluid 
stiffness follows a linear temperature dependence at low 
temperatures p s (T) ~ p s (0) — BT due to thermal exci- 
tation of quasiparticles at the Dirac nodes. The slope B 
is inversely proportional to the magnitude of the d-wavc 
gap. This expression implies a crude estimate of T c in 
a two dimensional system T c ~ p s (0)/(B + 2/tt), which 
is the point where the criterion for a Kosterlitz-Thouless 
transition is satisfied (p s (T c ) = 2T c /ir). Crucially, the 
transition temperature depends on the zero temperature 
stiffness, but also on the quasiparticle gap via the slope 
B. 

Coupling a layer of d-wave superconductor to a normal 
layer produces a d-wave proximity gap, which protects a 
super-flow of electrons in the normal layer. Fig. [1] is an 
illustration of the resulting low energy spectrum of the 
bilayer system. While the carrier density of the under- 
doped layer is small, proportional to hole doping of the 
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Mott insulator, that of the normal layer is much larger 
and is of the order of the total electron density. Conse- 
quently, the zero temperature superfluid stiffness is huge, 
consisting of the contributions from the two layers and 
is clearly dominated by the carrier density of the nor- 
mal layer p B (T = 0) « n s \jm\ + n s %Jm\. On the other 
hand the reduction of the stiffness with temperature is 
now much steeper than for a single layer. This is be- 
cause it is dominated by thermal excitation of quasipar- 
ticles at the nodes of the d-wave proximity gap, which 
is smaller than the pairing gap in the underdoped layer. 
At weak coupling between the layers, the proximity gap 
is given by A£ rox w Ak(?_i/£'i) 2 , where t± is an effec- 
tive inter-layer tunneling and E\ is a larger energy scale 
determined by the mismatch of the two Fermi surfaces. 
The resulting linear slope of the stiffness with temper- 
ature, is given to leading order in the small parameter 
tx/Ei by B 2 « (21n2/Tr)(v F2 /v A ){E 1 /t ± ) 2 - Here v F2 
is the Fermi velocity of the normal layer and «a the slope 
of the pairing gap of the underdoped layer at the node. 
This should be compared with the smaller slope in a pure 
underdoped material B\ = a 2 (21n2/7r) (vfi/va), where 
a is the effective electric charge carried by a current of 
quasiparticles. Thus the question of T c enhancement in 
the bilayer hinges on the competition between the in- 
creased zero temperature stiffness and reduced quasipar- 
ticle gap compared to the pure underdoped material. 

To understand the full dependence of T c enhancement 
on doping level and bilayer coupling one must go be- 
yond these leading order estimates. To do this we draw 
on the basic framework of slave boson mean field the- 
ory (SBMFT)ii^^. We also derive an effective semi- 
phenomenological theory of the bilayer, along the lines of 
Ref. [H, which keeps the spirit of SBMFT while avoid- 
ing some of its peculiarities. The resulting phase diagram 
shows significant enhancement of T c for moderate inter- 
layer tunneling over a wide range of doping levels. As 
a function of inter-layer tunneling, T c increases at first, 
but reaches a maximum at an optimal value of t±. The 
main factor in setting the optimal coupling is the anti- 
proximity effect on the pairing gap of the underdoped 
layer. 

We find another, somewhat more subtle, anti- 
proximity effect, which affects the zero temperature su- 
perfluid stiffness. In addition to the usual diamagnctic 
response, there is a paramagnetic correction at zero tem- 
perature due to mixing of quasiparticle wave-functions 
between the two layers. Accordingly, the zero tempera- 
ture superfluid stiffness is smaller than the independent 
contributions of the two layers, n s i/m\ + li^/mj, by a 
term proportional to (Ji — J2) 2 , where J/ = otivpi is 
the current carried by a quasiparticle located on layer 
I. We argue that measurements of this effect can lend 
insights into a long standing problem concerning the so- 
called quasiparticle charg e 12 ! 17 , or current carried by a 
quasiparticle in the cuprates. 

The rest of this paper is organized as follows. In sec- 
tion HD we define the microscopic model which serves as 



the basis for our theoretical analysis. In section IIIII we 
use self consistent slave boson mean field theory to obtain 
the temperature dependent phase stiffness and a phase 
diagram of the bilayer model. In section IIVI we derive a 
low energy effective theory of the response to an external 
vector potential starting the from the slave boson formu- 
lation. We then generalize the low energy theory to in- 
clude rcnormalizcd parameters for the zero temperature 
superfluid stiffness and the effective quasiparticle charge. 
In section HV CI we use the semi-phenomenological theory 
to predict T c enhancement in a putative LSCO bilayer 
composite system. In section IIVDI we use the effective 
theory to derive the paramagnetic correction to the zero 
temperature stiffness. Finally in section[V]wc summarize 
our main conclusions and discuss possible implications to 
recent experiments. 

II. THE MODEL 

Our starting point for theoretical investigation is the 
following model of a bilayer system: 

H = H 1 +H 2 + H ± 

Hi = -h J2 P Y\ a c 3 e + h.c.}P - (co + M) J2 c l c - 

(ij)a ia 

+J^[S i -S j --n i n j ] + ... 

m 

H 2 = -t a [dldja + h.c]- 

(ij)-a i 

H± = -t± 534r*o- + h.c}. (1) 

iu 

Here, the normal (highly overdoped) layer is modeled by 
the Hamiltonian H 2 of non interacting Fermions dl a on 
the square lattice. The underdoped layer on the other 
hand is modeled by the effective t — J Hamiltonian H\ , 
which takes into account the proximity of the Mott in- 
sulating state. P is the projection on the low energy 
subspace with no doubly occupied sites, and the dots 
represent possible additional terms. The energy offset 
eo is the single layer chemical potential that would set 
the correct hole doping of the underdoped layer in ab- 
sence of inter-layer coupling H±. When the two layers 
are coupled by H±_ they of course must share a common 
chemical potential /j,, which in general leads to charge 
redistribution between the layers. 

III. MEAN FIELD PHASE DIAGRAM 

In this section we obtain the phase diagram of the 
bilayer model using the slave boson mean field theory 
(SBMFT^Liii&Ii. This is the simplest theory, that gives 
a BCS like superconductor with a large Fermi surface 
but with low superfluid density, which scales as the hole 
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doping. Thus for a single layer T c is controlled by the 
zero temperature supcrfluid stiffness, rather than by the 
pairing gap. 



where 5 = \j2x/(l + x). The quadratic bilayer (varia- 
tional) Hamiltonian in momentum space is then 



H M f = £[£ik + 6k + *£/ik*k], 



(5) 



A. Slave boson mean field theory for a bilayer 

Before presenting the bilayer calculation, let us briefly 
review the standard slave boson approach for a single 
underdoped layer. The electron creation operator is 
represented as a composite of a fermionic spinon and 
a bosonic holon operator c\ a = &i/ i(T . The redun- 
dancy of representation is removed by the local constraint 
bl^i +^2 a f tafia = 1, which can be implemented exactly 
by a U(l) gauge field. The core approximation of the 
mean field solution is that at least in the superconducting 
phase, both the holon and the gauge field are condensed. 
This allows to replace the operator bi by the number 
[2x/(l + a:)] 1 / 2 = y/S and implement the constraint only 
on the average^. At this stage the approximate Hamil- 
tonian H'(x) is written in terms of the fermion spinon 
operators only, and acts in an unrestricted Hilbert space, 
however it is still quartic. The second approximation 
consists of a standard mean field solution of H', whereby 
one seeks the best quadratic approximation to it of the 
form 



H = -h5j2(j eAii fIfj + h 



<</} 



c. - e 



■ £(A k / T yi,-k^) - E XkfLf*, (2) 



k<7 



where Ak = A(cosfca; — cosfe^) and Xk = xi cos ^x + 
cos k y ). Here we introduced a coupling to an external 
vector potential through the phases eAy , which will later 
facilitate calculation of the supcrfluid density. Note that 
the electromagnetic vector potential couples only to the 
charged holon field. Condensation of the holon leads to 
effective coupling to the fermion field in the kinetic en- 
ergy term. The parameters A and x & r e determined using 
a general thermodynamic variational principle by mini- 
mization of 



(H' - H, 



0/0) 



(3) 



Fq is the free energy implied by the trial Hamiltonian 
H , and () denotes a thermal average generated by H - 
The chemical potential eo is determined by resolution of 
the average constraint equation: (foifd) + % = 1. 

We now move on to include the inter-layer coupling. At 
this point charge can be redistributed between the layers, 
changing the doping levels of the two layers from x and 
y in absence of the coupling to x and y. The quadratic 
inter-layer tunneling Hamiltonian is given by 



Hxo = -t x VtJ2lfl 



(4) 



where ^ 



7kT'/-kl! rf kT' d - k i) and 



/6k A k i x ^ 

A k -£i,-k -to. 

t± o 6k o 

V o -i± o -6,-k/ 



(G) 



with tj_ = t± Vs. In the absence of external fields 

£ik = — (2Sti + x) (cos k x + cos k y ) — /j. 
6k = -2£ 2 (cos k x + cos k y ) - /i + e 
Ak = A(cos k x — cos k y ). 



(7) 



The parameters A and \ can be determined again by 
solving the variational equations, supplemented by the 
two number equations for the additional unknowns x and 
/x: 



(n d ) 



x 

y 



1 — y + (x — x). 



(8) 



Clearly, a proximity gap will be induced in the normal 
layer due to the coupling with the underdoped supercon- 
ducting layer. The charge carriers in the second layer 
will then contribute to the supcrfluid density. 



B. Superfluid density and T c 

We shall obtain the critical temperature of the bilayer 
by computing the temperature dependent supcrfluid stiff- 



Ps(T) 



1 d 2 F 



ndA 2 



(9) 



A=0 



Here F is the free energy, A is an externally applied trans- 
verse vector potential and VL is the volume of the system. 
The critical temperature is then given by the condition 
for a Kosterlitz-Thouless transition p s {T c ) = (2/-k)T c . 

The dominant contribution to the reduction of super- 
fluid stiffness at low temperatures in a d-wave supercon- 
ductor is the paramagnetic response due to thermally 
excited quasiparticles in the gap nodesii. This leads to 
a linear decrease of p s with T. Because the paramag- 
netic response amounts to a current-current correlator, 
the slope dp s /dT is proportional to the square of the ef- 
fective electric charge carried by a current of quasiparti- 
cles. Here we encounter a possible pitfall of the mean field 
theory. Within plain SBMFT, the quasiparticlc charge 
is proportional to the doping x. However, experiments 
seem to point to a fairly doping independent value of this 



4 



1 

0.8 
| o 0.6 
H ° 0.4 

0.2 


1.5 



(a) , " " 
s 

/ 

/ 


N 

S 

S 

t ± =0.5t 


/ 

/ 


— ► =0.3t - 


/ f 

/ y 


t=0.1t 

_L 







0.05 0.1 0.15 0.2 0.25 0.3 0.35 

x 



0.5 



(b) ^ 






t =0.5t / / 


>5 




t,=0.3t / / 

A 






t ± =0.1t /' 






/ 









-0.05 



0.05 



0.1 

x 



0.15 



0.2 



0.25 



FIG. 2: Phase diagram of the bilayer system from microscopic 
theory. The critical temperature T c (normalized by its maxi- 
mal value for two identical underdoped layers, T™q*) vs. dop- 
ing of the underdoped layer. The doping level of the metallic 
layer is y = 0.35. The dashed line is the result for two iden- 
tical underdoped layers, (a) "Bare" SBMFT calculation. No 
enhancement of T c . (b) SBMFT with renormalized quasipar- 
ticle charge of a± — 0.5 in the underdoped layer. Maximal T c 
enhanced by ~ 40% in the heterostructure. Optimal doping 
shifted down, consistent with experiment of Yuli et. al. jlOi ]. 



paramete r 18 ! 19 . The experimental results can be repro- 
duced by an effective theory, which maintains the spirit 
of SBMFT, but assigns a phenomenological value to the 
quasiparticle charge^. 

In our analysis of the two layer system we consider the 
bare SBMFT as well as a theory with a phenomenological 
quasiparticle charge renormalization. The main result 
of the mean field calculation is a phase diagram of the 
bilayer heterostructure. 

Fig. [2fa) displays T c as a function of the doping x 
computed using the bare slave boson theory for various 
values of the inter-layer tunneling t±. No enhancement 
of T c relative to a pair of identical layers is found for the 
relevant range of parameters 0.2 < J < 0.7 and < t± < 
1. 

In marked contrast, we do find a significantly enhanced 
T c in a modified SBMFT which allows for a phenomeno- 



logical quasiparticle charge renormalization independent 
of the doping. To compute the second order response © 
to an external vector potential A, the mean field Hamil- 
tonian © is expanded to second order in Aij 



H(A) = H(0) 



(ij) 



A 2 



(10) 



with the paramagnetic current operator jij — 
ttiSeJ^Miofj" ~ fjafi") and tne average kinetic en- 
ergy per bond k i3 = -Me 2 Y,Mlfj° + f] a fia)- In the 
plain SBMFT approach the charge in the current opera- 
tor jij is renormalized by a factor 5 cx x. It is this renor- 
malization that leads to a strong doping dependence of 
the slope dp s /d,T, at low temperature*^, which disagrees 
with experiment a 18 i 19 . As a possible cure of this artifact 
within the microscopic theory, we replace the factor 5 in 
jij by a doping independent number a, which in principle 
should be determined experimentally. This is equivalent 
to introducing an effective quasiparticle charge of mag- 
nitude ae to all physical properties involving the quasi- 
particle current, as suggested in Refs. [ill and[H|. In our 
model we apply this renormalization only to the under- 
doped layer since the normal layer is approximated sim- 
ply by non interacting fcrmions. We choose a value of 
a sa 0.5, which reproduces the dome-shaped T c (x) phase 
diagram with T c jt\ having the right order of magnitude. 
For this renormalized value we obtain enhancement of 
T c of the bilayer compared to two identical underdoped 
layers with the same a. 

The result of this calculation is presented in Fig. &b)- 
The optimal doping level of the underdoped layer of the 
heterostructure is seen to be around x = 0.1, which is 
well below the optimal doping of the single layer and 
consistent with the result of recent experiments^. 

The enhancement of T c stems from the combination of 
a large carrier number donated by the normal layer and a 
large pairing gap induced by the proximate underdoped 
layer. Because the proximity gap is smaller than the 
original pairing gap, the reduction of stiffness with tem- 
perature is also steeper in the heterostructure. However, 
for the modified SBMFT with renormalized quasiparticle 
charge, the increase in slope dp s /dT is not large enough 
to offset the enhanced zero temperature superfluid stiff- 
ness. By comparison, in the plain SBMFT the slope of 
stiffness versus temperature for the pure underdoped ma- 
terial is much smaller because of the small quasiparticle 
charge. The increase in slope dp s /dT on going from a 
pure underdoped material to a heterostructure is con- 
comitantly more extreme. For this reason we see an 
enhancement of T c only for SBMFT with renormalized 
quasiparticle charge. 

We note that the computation of superfluid stiffness 
© with both plain and modified SBMFT is carried out 
fully self consistently. Thus it captures non linear contri- 
butions to the temperature dependence of the superfluid 
stiffness (but of-course not contributions from phase fluc- 
tuations). 
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The degree of T c enhancement as a function of the 
inter-layer tunneling is plotted in Fig. [3] We note that 
the minimal inter-layer tunneling required to obtain such 
an enhancement (t±_ ~ i/4) appears rather too large to 
serve as straightforward model of the bilayer experiment 
of Yuli et. al [lfj. This issue will be discussed further 
in section IIV CI Moving to still larger t± we observe a 
maximal enhancement of T c at t± ~ 0.5i for which the 
enhancement may be as large as 40%. The optimum 
value of t± occurs where the proximity gap becomes of 
order of the superconducting gap (see inset of Fig. [3]) . At 
this point the superconducting gap cannot increase any 
further and the anti-proximity effect of the normal layer 
on the superconducting one takes over. 

Before closing this section we point out another inter- 
esting effect in the response of the bilayer heterostructure 
as compared to the single layer (or a pair of identical 
layers). In the single layer mean field theory the para- 
magnetic response appears only due to quasiparticles at 
finite temperatures, whereas the zero temperature super- 
fluid stiffness is independent of the effective quasiparticle 
charge. It is not so in the coupled bilayer heterostruc- 
ture. In this case a quasiparticle carries a different cur- 
rent depending on whether it resides on the top or bottom 
layer. This allows redistribution of quasiparticles in the 
ground state in the presence of current and leads to a zero 
temperature paramagnetic contribution to the supcrfluid 
stiffness. 



We will show in section HVDi that the reduction 
of the stiffness at zero temperature is proportional to 
(aiV.Fi — CY2VF2) 5 where ai and vfi are the effective 
quasiparticle charge and the Fermi velocity in layer /. 
Note that aiVpi = J; is the current carried by a quasi- 
particle in layer I. Fig. [5] shows the zero temperature 
stiffness of the bilayer within the self consistent mean 
field calculation as a function of the quasiparticle charge 
renormalization. We see indeed that the correction is 
negative, and quadratic in J\ — Ji- 
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FIG. 3: Optimal inter-layer tunneling. Plotted is the critical 
temperature T c vs. the inter-layer coupling t±/t (t = ti = 
£2). Different curves correspond to different doping levels of 
the underdoped layer and all curves are normalized by the 
critical temperature of two underdoped layers of the same 
doping level, T c ,o(x). The doping of the metallic layer is y = 
0.35. Inset: The self-consistent gap A and the proximity 
gap A prox as calculated from the bilayer energy spectrum, 
plotted vs. the inter-layer tunneling t±/t for doping x — 0.11. 
This suggests that the optimal value of t± is determined by 
the point where the anti-proximity effect on the gap of the 
underdoped layer overtakes the proximity gap in the normal 
layer. 



A. Quasi-particle spectrum 

In the absence of inter-layer coupling, the lower diag- 
onal block in the microscopic Hamiltonian ^ describes 
gapless particle and hole excitations near the Fermi sur- 
face of the metallic layer — 0). The primary effect 
of the coupling is to open a proximity gap in the metallic 
layer. This is captured nicely by the low energy effective 
Hamiltonian derived by second order degenerate pertur- 
bation theory. To this end it is convenient to rewrite the 
Hamiltonian (El) in terms of its 2X2 blocks 



IV. EFFECTIVE LOW ENERGY THEORY 



To clarify the mechanisms of T c enhancement and fa- 
cilitate generalizations that are less dependent on a par- 
ticular microscopic model it is worthwhile to derive a 
low energy effective theory for the bilayer system. This 
will be done perturbatively in the inter-layer coupling in 
subsection IIV Al We shall also derive the effective cou- 
pling of the external field to the low energy Hamiltonian 
for the sake of computing the superfluid stiffness. The 
low energy theory, with parameters extracted from bulk 
samples, will then be used to construct a phase diagram 
of the bilayer system for given values of the inter-layer 
coupling. 



hir = 



Mk) 



V 



(11) 



and treat V as a perturbation. We note that since the 
effective inter-layer tunneling t±_ is proportional to 
at sufficiently low doping levels a perturbative treatment 
may be justified even if the bare tunneling tj_ is not very 
small. 

The low energy physics is dominated by excitations 
near the Fermi surface of the metallic layer (layer 2) and 
near the nodal points of the superconducting layer (layer 
1). The effective Hamiltonian for the metallic layer near 



its Fermi surface (i.e. 
standard way2£ 



t(2) 



( 2") 

— £^ ) is obtained in a 
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FIG. 4: Paramagnetic reduction of the zero temperature stiff- 
ness. p s (T — 0) of the bilayer system (doping levels x = 0.11 
and y = 0.35) is plotted against (Ji — J2) 2 /t 2 , where t = ti = 
i2 is the bare in-layer tunneling and J; = ctiVfi are the quasi- 
particle currents of the two layers. Here vfi are the Fermi 
velocities at the nodes in the respective layers. We used the 
quasiparticle charge a.2 = 1 for the normal layer and varied 
u\ of the underdoped layer. 



on the underdoped film, for which no enhancement of T c 
was found in Reffiol. 

In addition to (fT3"j). the complete low energy Hamilto- 
nian also includes the Dirac quasiparticles of the origi- 
nal superconducting layer, which are now slightly renor- 
malized by degenerate perturbation theory (near i?i k ~ 
--Elk), 

h eS _ ( £l,k — ij_/£2k A k \ 

^"V A k -6 ) -k + il/6ky'- 

The low energy effective theory captures correctly 
properties related to quasiparticle excitations at low tem- 
perature. To get a full picture of zero temperature prop- 
erties we should include perturbative corrections to all 
negative energy states, including those far below the un- 
derlying Fermi surfaces of the two layers. Such correc- 
tions will be discussed in section TlVDi where we analyze 
a unique paramagnetic contribution to the stiffness at 
zero temperature. 

B. Response to transverse vector potential 



Note that the effective Hamiltonian is energy dependent, 
and therefore not really a Hamiltonian. This is because 
it is defined through the resolvent operator projected to 
the lower right block 

G 22 (E) = P 2 [E- h]- x P 2 = [E- hf{E)]~\ (12) 

The energy dependence will be important below when we 
consider response to external fields. However, for now, 
since we are only interested in the low energy spectrum, 
compared to the separation between blocks, we may ne- 
glect the energy dependence to leading order in degener- 
ate perturbation theory and obtain 

h cs /&,k-(*J-/£ik) 2 £i,k («i/£ik) 2 A k \ 
2 ' k V (W£ik) 2 Ak -6,-k+(i±/£ik) 2 £i,- k ; ' 

(13) 

Here -Ei k = \/€ 2 k + A| is the energy of a quasiparti- 
cle of the superconducting layer at a wave vector k near 
the Fermi surface of the normal layer. We see that a 
small t± leads to a proximity gap A£ rox = (t±/En ti ) 2 A] <: . 
The proximity gap inherits the d-wave symmetry from 
the pairing gap A k of the superconducting layer, but 
is suppressed in magnitude. An important observation 
is that the degree to which the proximity gap is sup- 
pressed is highly sensitive to the Fermi surface matching 
between the two layers. For highly matched Fermi sur- 
faces the energy denominator i?i k is very small compared 
to the full bandwidth 4t±. In this way it is possible to 
gain significant enhancement in T c with relatively small 
inter-layer tunneling t±. It is interesting to note that 
between a highly overdoped layer with hole concentra- 
tion x ~ 0.35 and an underdoped layer, as seen in angle 
resolved photocmission (ARPES) experiments, is impres- 
sively good^i. This is not so in the case of gold deposited 



To compute the superfluid stiffness within the low en- 
ergy effective theory using formula @ , we need to derive 
the renormalized coupling to an external vector potential 
in the effective Hamiltonian. This is accomplished by car- 
rying out the renormalization scheme outlined above in 
the presence of a field, while keeping terms up to sec- 
ond order in A throughout. For simplicity let us take 
A = Ak. 

The coupling of the microscopic Hamiltonian © to the 
vector potential, up to second order in A, is given by 

h t (k, A) = h, (k, 0) - J, • A + 1 K X (k) A 2 (14) 

where the index I = 1,2 refers to the two layers, and 
A';(k) = d 2 fiCT3 = £;"c3 are the kinetic energy operators 
due to motion along the axis defined by A (in our case the 
x axis). 3 1 = aiVpicro, is the electric current operator on 
the layer I, with ai the quasiparticle charge on that layer 
and vpi the Fermi velocity. The off diagonal block V of 
the Hamiltonian does not couple to the electromagnetic 
field. 

Now following the same steps as above we can elimi- 
nate the coupling between the blocks and obtain an ef- 
fective Hamiltonian valid near the Fermi surface of the 
metallic layer 

hf{E,A) = h 2 -3 2 A + V^[E-h 1 +J 1 -A--K 1 A 2 }- 1 V 7 

(15) 

where we have dropped the argument k for notational 
simplicity. We note that the energy E should be under- 
stood as a solution to the equation det(E — h2 S {E)) = 
for the poles of (fT"2"|) . Therefore the energies implicitly de- 
pend on the external field A. To zeroth order in t±/Ei 
we have E(A) = E(0) - J 2 • A + isgn(£(0))£ 2 ' A 2 . We 
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must keep the A dependence since we are interested in E\ in this regime. In this way we obtain the effective 
the response to the external field. However we may still Hamiltonian 
neglect the constant E(0), which is much smaller than 



hf (A) = h 2 - J 2 • A + h< 2 A 2 - V*[hi - (Ji - J 2 ) ■ A + i(sgn(E) $a - ga 3 )A 2 ]- l V 

= h 2 - V^hfV - J 2 • A + X -K 2 A 2 - (Ji - J 2 ) • A 

I 



^ ) ( hf(J 1 J 2 y - (sgn(E)e> - £V 3 ) ) ^. 



(16) 



This is still energy dependent because of the term sgn(E). 
However, of the two terms quadratic in A in the last line, 
the first is larger by a factor ~ At\/E\ ~ it /8k (see ap- 
pendix for detailed explanation), where 5k is the mis- 
match between the Fermi surfaces of the two layers at 
the nodes. In other words the first term is strongly en- 
hanced by good Fermi surface matching, which is indeed 
observed by in ARPES experiments done with samples 
of varying doping levels^. Specifically, for a bilayer with 
underdoped layer at x = 0.07 — 0.15 and overdoped layer 
doping y — 0.35 the estimated ratio is ti/Ei > 5. We 
therefore neglect the energy dependent term and obtain 
an effective Hamiltonian 



hf{A) = h c 2 \-3 



eff 



reff 



A + -KfA 2 
2 2 



(17) 



with 



Jf = J2 + ( (Ji - Ja) 

f'tV^^(Ji-J 2 ) 2 .(18) 



£1 



Ei 



The effective Hamiltonian hf 1 (A) valid near the Dirac 
nodes of the underdoped layer is derived in the same way, 



hf (A) = hl% - Jf • A + -KfA 2 (19) 



where 



of the stiffness with the temperature. Because this 
reduction is induced by thermal excitation of low energy 
quasiparticles at the gap nodes, it can be computed 
using the effective low energy theory. This will be done 
in the next subsection. We note that there are also zero 
temperature corrections to the stiffness due to coupling 
between the layers. These are somewhat more subtle 
and will be considered in subsection IIVD1 



C. Temperature dependent phase stiffness 

In the low energy effective Hamiltonian we have 
achieved effective decoupling of the two layers. There- 
fore, the contributions of the stiffness due to each layer, 
within this theory, can be added separately and they 
must each be non negative 



PsiT) = max(pi,0) + max(/) 2 , 0). 



(21) 



At zeroth order in the inter-layer coupling, the zero 
temperature stiffness is simply the sum of the contribu- 
tions of the independent layers, pi(0) = n s i/m^. The 
leading temperature dependence of pi is a linear reduc- 
tion in temperature due to a paramagnetic contribution 
from thermally excited quasiparticles in the nodes of Kf 
and h 2 s . This contribution can be calculated exactly as 
in Ref. [ll[ using the effective Dirac Hamiltonians hf 
and the respective quasiparticle currents jf : 



Jf = J 1 



(Ja- Ji) 



S PltPm = - f £ (Jf)" n F (k) (1 - n F (k)) . (22) 



Kf = A' 1 -2-^ ( t 3 (Ji-J 2 ) 2 . 

S2 



(20) 



The phase stiffness we wish to compute can be divided 
into two parts. First is the zero temperature superfluid 
stiffness, which to leading order in the inter-layer 
coupling is given simply by the sum of contributions 
from the two layers. Second, is the linear reduction 



We carry out the integration using the density of states 
of the respective layers, 



ME) 
ME) 



E 



2itvfiVa 
E 

2ttvf 2 va 



(23) 
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where va = va^/E 2 and vf2 = vf2 — VFit\/ 'Ef, (fT8|) and ([20)) we obtain the low temperature contribu- 
the proximity induced gap and Fermi velocities near the tions to the phase stiffness due to each of the layers 
Fermi surface of the metallic layer. Taking Jf ff from Eqs. 



Pi 



Pi 



2 In 2 2 v F1 Mn fP ± 
pi(0)-T a t +0 -± 

7T V A V^2 



21n2 



I' A 



E 



VF1 ! 2 , o 

(a 2 + 2aia.2 

VA " 



(24) 



Clearly the dominant term in the temperature depen- 
dence of p s , is due to thermal excitation of quasiparticles 
in the proximity-induced Dirac cones of the metal (layer 
2). This term scales as vf2/va = {vF2/vA)(Ei/t 2 L ). 
Terms of order 1 in the dimensionless inter-layer tun- 
neling t±/Ei are due to excitations in the original Dirac 
cones and to quasiparticle mixing between the layers. 

Our next step is to estimate^ T c of the bilayer using 
the formulae (f2"Tj) and (|2~i|) . These formulae are expressed 
mostly in terms of phenomenological parameters, which 
may in principle be extracted from experiments with bulk 
samples. 

We estimated the needed parameters using the fol- 
lowing information: (i) The zero temperature stiffness 
of the underdoped layer pi(0) was taken from data in- 
terpolation of penetration depth measurements— (ii) 
In the underdoped regime dpi/dT r* — 1 almost inde- 
pendent of doping— (iii) va is determined from the 
maximal gap extracted from the leading edge shift in 
ARPES 2 ^ by assuming a pure d-wave gap function A k = 
A(cos k x — cos k y ). (iv) The Fermi velocity of the under- 
doped material is taken from ARPES measurement s 21 ' 26 , 
which give ~ 1.8 eV — A almost independent of the dop- 
ing within the underdoped regime, (v) The effective 
quasiparticle charge in the underdoped regime, ot\ (x) = 
[{Tr/21n2)(v A /v F i)dpi/dT} 1 / 2 , is then fully determined 
by (ii)-(iv). (vi) The zero temperature stiffness of the 
metallic layer is estimated as the average kinetic energy 
per bond, P2(0) = £2(1 — y), where we take y = 0.35. 
The hopping £2 ~ 300 meV is taken from the band 
structure determined by ARPES measurements of LSCO 
samples— The ratio of the two Fermi velocities is seen 
to be vf2/vfi ~ 1-5. 

The one parameter that cannot be extracted from 
such experiments is the dimensionless inter-layer cou- 
pling t±/E\. We remind the reader that -Ea(k) is the 
energy of a quasiparticle of the superconducting layer at 
the wave- vector k near the Fermi surface of the metallic 
layer. It therefore depends crucially on the distance 8k 
between the two Fermi surfaces and can be approximated 
as Ei ks Skvpi- In principle 8k may be extracted from 
ARPES experiments, such as Ref.[2l|]. For underdoped 
LSCO layer with hole concentration between x = 0.07 
to x = 0.15, matched with a highly overdoped layer 



x = 0.35, we extract E\ ss 60 meV <C t. It is there- 
fore possible, in principle, to obtain a sizable proximity 
effect, even for inter-layer tunneling substantially smaller 
than t as long as t± is not much smaller than E\. 

Fig. [5] shows the phase diagram of LSCO bilayers es- 
timated using the phenomenological theory described in 
this section. In the underdoped side T c is controlled by 
the temperature dependence of the supcrfluid stiffness of 
the bilayer as given by Eqs. (|2"4")l and (f2"Tj) . The value 
of T c is determined by the criterion p s (T c ) = (2/n)T c for 
a Kosterlitz-Thouless transition in the two dimensional 
interface layer. The result for T c (x) of the double layer 
systems for two values of the inter layer coupling is given 
by the solid lines in the figure 2 ^. These lines are cut off 
by the dashed curve which is the estimate of T c more 
appropriate to the overdoped side of the phase diagram. 
There, the pairing gap becomes smaller than the energy 
scale set by the superfluid stiffness and therefore the gap 
sets the scale for T c , which can be estimated by the BCS 
relation Tc ~ A(x)/2. The value of the gap as a func- 
tion of doping is taken from an interpolation of ARPES 
data—. For comparison we also show the transition tem- 
perature measured in bulk LSCOSiL (black circles). 

For inter-layer tunneling t± > t/5 we see a significant 
enhancement of T c compared to the bulk transition tem- 
perature. This is one of our main results. Furthermore 
the optimal doping level is shifted down compared to 
the bulk optimal doping, in qualitative agreement with 
experiment—. We point out that for inter-layer tunneling 
t± = t/5 the perturbative parameter is (t±/Ei) 2 ~ 0.4, 
justifying the expansion in (|24[) . However, we note again 
that the inter-layer tunneling required to achieve the en- 
hancement of T c is rather large to directly explain this 
experimen t 28 ' 29 . 

One possible explanation of the experimental result, 
in-line with our analysis, is that in reality the interface 
layers share dopants, such that each layer is an inhomo- 
geneous mixture of underdoped and overdoped puddles. 
This is a natural scenario in the samples of Ref. ITol . in 
which the interface is not atomically sharp and it was 
shown to consist of facets of the two material compo- 
nents. However inhomogencous doping of the interface is 
plausible even in the atomically sharp interfaces of Ref. 
9. Indeed these authors mapped the doped hole distri- 
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0.25 



FIG. 5: Phase diagram of bilayer LSCO from the phenomeno- 
logical theory. Plots of the critical temperature for a bilayer 
with different inter-layer tunneling t±/t computed using the 
phenomenological approach of section HV CI (solid red curves). 
This is compared to the measured critical temperature in bulk 
LSCO (black circles) taken from Ref. [24| . The dashed line is 
a linear interpolation of data for A(x)/2 (anti-nodal gap)2^, 
taken as an estimate for the mean-field critical temperature. 



bution along the c-axis using resonant X-ray scattering 
and found that the interface doping is approximately the 
average of the nominal doping levels of the two material 
components 2 ^. This was explained by a simple theory of 
electrostatic screening. The hole distribution within the 
interface plane was not mapped, but it is highly likely to 
be inhomogeneous given the random dopant distribution. 

The essential competition of effects that determine T c 
in a inhomogeneous layer is expected to be the same as 
discussed above. The large superfluid density donated 
by the proximity gapped overdoped regions counters the 
steep reduction of stiffness with temperature due to the 
smallncss of the proximity gap. Most importantly, now 
that the two 'phases' are intertwined in the same layer, 
the effective coupling between them can be much larger. 
Our analysis of two homogeneous layers with substantial 
coupling between them could then be viewed as a crude 
effective description of the inhomogeneous system. 



D. Zero temperature paramagnetic response 



(fP7|) and (fT9]) . In the effective hamiltonian of each of 
the layers the quadratic coupling to a vector potential 
is renormalizcd down, at second order in the inter-layer 
tunneling, by a factor proportional to (a±vi — a^Vi) 2 . 
This is not the full contribution to the zero temperature 
stiffness, which consists of the response of the ground 
state energy to the vector potential. The ground state 
energy, in turn, involves a sum of all the negative energy 
solutions of ([6]). It is therefore not enough to compute the 
contribution from the low energy excitations, encoded by 
the effective Hamiltonians (fTT)| and (jTTJ). As opposed to 
the calculation of the temperature dependence presented 
above, here we must also account for the contribution of 
the negative energy solution of the high energy excitation 
branch. That is, the energy E\-(k.,A) at wave-vectors 
near the Fermi surface of the metallic layer (layer 2) and 
E2- (k, A) near the nodes of the undcrdoped layer (layer 
I)- 

Thus, the contributions to the superfluid stiffness from 
the wave- vectors near the fermi surfaces of the two layers 
are given by: 



Pa(0) 



kg{a} 



dA 2 



eff\ 



(25) 



Here the two layers are denoted by a = 1,2, while a = 
2, 1 denotes the other layer. K^f are given in (jT5J) and 



(|20)) and (...)o denotes a ground state expectation value. 
Using the same procedure as outlined in section lTVB| but 
applied to the large negative energy solutions, we get 



d 2 E a _(A) 
dA 2 



set 



E 3 



(Ji 



J 2 ) 2 . 



(26) 



Finally, using (|25p and (|26|) and interpolating to all wave- 
vectors we obtain the correction to the zero temperature 
stiffness to order (t±/E a ) 2 



Sp(0) ~ —(aiVpi - a 2 VF2) 



2t 2 j_ sin 2 k x 



2 ^(£ik + |£2k|) 3 



(27) 



2k 



-Elk 1 6k 



So far we were concerned with the variation of the stiff- 
ness with temperature. In this section we point out an 
interesting zero temperature effect of the bilayer coupling 
on the superfluid stiffness. Specifically, the superfluid 
stiffness of the coupled double layer system is smaller 
than the summed stiffness of the individual layers. 

One way to see this effect is by inspection of the ef- 
fective hamiltonian of the double layer, as given by Eqs. 



This is added of course to the zeroth order stiffness of 
the two layers ~ n s \jm\ + n S 2jm\ (see Eq. IA4|) . 

To gain better understanding of the zero temperature 
paramagnetic correction and the processes involved, it is 
worthwhile to derive it from a diagrammatic approach. 
At second order in t±, the diamagnetic and paramag- 
netic corrections to the superfluid density are given by 
the following diagrams 
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Spdi 



Sp pa 




+ 2 x 




(28) 



r 



Here the labels 1 and 2 denote the bare Nambu Green's 
functions of the isolated layers 1 and 2 respectively and 
a vertex with a dashed line denotes an inter-layer tun- 
neling process. In adding up the contributions of the 
paramagnetic bubble diagrams above, we directly obtain 
the result (|2"T|) . Note that the first two diagrams describe 
renormalization of the Green's function of each layer due 
to virtual hopping of electrons out of it into the other 
layer. The third diagram is a vertex correction describ- 
ing indirect coupling to the vector potential via hopping 
to the other layer. The correction to the stiffness of the 
same order in ij_ coming from the diamagnetic diagrams 
is suppressed by a factor of order Sk/kp ■ It can therefore 
be neglected in the case of good Fermi surface matching. 

In the appendix we present a more complete derivation 
of the result (|27|) using degenerate perturbation theory. 
In addition we show there that the paramagnetic cor- 
rection behaves as (aivi — a2V2) 2 at all orders in t±. 
The same conclusion also emerges from the results of the 
self consistent mean field calculation (described in sec- 
tion [TTTJ) . Fig. 2] shows a quadratic dependence of the 
zero temperature stiffness on ( J\ — J 2 ) when the effective 
quasiparticle charge of the undcrdoped layer is varied. It 
is interesting to note the weak dependence of p s (0) on 
t± at the point Ji — 3 2 = 0. This arises from the small 
diamagnetic term, which we neglected in the analytic cal- 
culation. 

The paramagnetic correction at zero temperature can 
be used to measure the current carried by a quasiparticle, 
and specifically its doping dependence. As mentioned 
above, this property also comes up in the temperature 
dependent stiffness, and its dependence on hole doping 
has posed a long standing puzzle (see for example Ref. 
[ItJ )■ Since these measures of the quasiparticle charge 
(or current), are model dependent it is useful to have an 
independent probe, such as the zero temperature param- 
agnetic effect in a bilayer. 

To define a concrete experiment along these lines it is 
simpler to consider a bilayer or heterostructure consisting 
of two types of underdoped layers with a small mismatch 
in doping \x — y\ <C x,y. The experiment involves com- 
parison between the supcrfluid stiffness measured for the 
heterostructures to that of the pure materials. Using the 
diagrammatic perturbation theory to second order in the 



inter-layer coupling we obtain the reduction of the zero 
temperature stiffness in the heterostructure: 



(29) 



8p{0) ~ -(aivpi - a 2 v F2 ) 2 V /^-L sm 



£lk$2k — AikA2k 
-Elk-E^k 



Note that here the electron operators of both layers are 
renormalizcd and therefore the effective interlayer cou- 
pling is t± = y // S(x)S(y)t±. The main contribution to 
the sum (|29[) is from wave- vectors between the underlying 
fermi surfaces of the two layers. Thus in the limit of good 
fermi surface matching between the layers (dk <C kp) we 
obtain 



5p(0) oc 



xy(ai - a 2 ) 



2 4-2 



5k 2 (x, y) 



(30) 



where we assumed that in the undcrdoped regime the 
fermi velocity is doping independent and thus vp\ = 
i'F2 = vf- The fermi wave vector depends very weakly 
on the doping and we denote by kF the average value of 
the two layers. In addition we plugged S(x) oc x which is 
valid at low doping levels. If we assume that 5k oc x — y 
(as should be expected from the Luttinger theorem), a 
measurement of the paramagnetic reduction at T = as 
function of both x and y can reveal the doping depen- 
dence of the effective quasiparticle charge a in the undcr- 
doped regime. Such a measurement will distinguish be- 
tween the following scenarios: (i) a(x) = a independent 
of doping, in which case there will be no paramagnetic 
reduction. Note that the small diamagnetic correction 
that will survive in this case is positive and therefore 
cannot be mistaken with the paramagnetic correction, 
(ii) a(x) ~ x as implied by slave boson mean field the- 
ory, which will result in a finite paramagnetic reduction 
that scales as x 2 to leading order in y — x. (iii) a{x) 
has some other doping dependence, leading to a more in- 
volved doping dependence of the paramagnetic response. 
In general, if a(x) depends on the doping as x 1 (7 ^ 0), 
then the leading doping dependence of the paramagnetic 
reduction is x 21 . 
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V. SUMMARY 

In this paper we showed that significant enhancement 
of T c in cuprate heterostructures is possible under real- 
istic conditions and provided a possible explanation for 
recent measurements on LSCO bilayers by Yuli et. al. 
[Toj | . Our analysis indicates that the conditions under 
which such enhancement of T c can occur depend cru- 
cially on the evolution of the superconducting state with 
underdoping on approaching the Mott insulator. In par- 
ticular, the effect is sensitive to the way in which the 
phase stiffness and the current carried by quasiparticles 
are renormalized as a function of the doping level. Such 
questions, pertaining to the fundamental nature of super- 
conductivity in the cuprates, are not yet fully resolved, 
and we proposed that further experiments with cuprate 
heterostructures can shed new light on these issues. 

The essential idea of T c enhancement in 
heterostructures^ is based on the observation that 
the pairing scale in the underdopcd cuprates is high, 
and T c is limited by the low supcrfluid density in these 
materials^. By inducing a proximity gap in a nearby 
metallic layer, the large density of charge carriers in 
that layer is harnessed to the total superfluid response. 
However because the proximity gap is typically much 
smaller than the original gap, the reduction of the 
superfluid density with temperature is also much steeper 
in the heterostructure. And so, the question whether 
T c can in fact be enhanced in this way is more delicate, 
and sensitive to the nature of superconductivity in the 
under doped material. 

To address this question we used a microscopic ap- 
proach based on the slave Boson mean field theory, as well 
as a semi-phcnomcnological theory of the doped Mott 
insulator. Straightforward slave boson mean field the- 
ory showed no enhancement of T c in the bilaycr. Inter- 
estingly however, this failure is directly tied to the well 
known shortcoming of the mean field theory in describ- 
ing the temperature dependent phase stiffness (see e.g. 
[171] ). One can generalize the low energy theory derived 
from the microscopic approach to include renormalized 
parameters for the zero temperature stiffness and the ef- 
fective charge of a quasiparticlc such that it reproduces 
the observed response in bulk samples^. Using such a 
phcnomenological theory for the bilayer we found that 
T c enhancement can be achieved for inter-layer tunnel- 
ing of order t/b or larger. This value is in excess of the 
bare inter-layer tunneling in LSCO bilayers, such as those 
investigated in Ref. [13] ■ We proposed that this discrep- 
ancy may be resolved if each of the layers at the interface 
is in fact an inhomogeneous mixture of undcrdoped and 
overdoped material (e.g. as a result of dopant migration). 
In this case our bilayer model with moderate coupling 
~ t/b can be viewed as a crude effective description 
of the inhomogeneous interface. 

We note that the analysis performed in this paper uses 
a completely homogeneous model. It docs not include for 
example stripe or density waves structures. The existence 



of such structures therefore does not appear to be crucial 
for obtaining an enhanced T c . From this point of view 
the fact that the maximal enhancement seen in Ref. fT3 | 
was close to 1/8 may be coincidental. 

Finally, we pointed out a unique paramagnetic contri- 
bution to the zero temperature phase stiffness of a bilaycr 
system. The paramagnetic reduction of the zero temper- 
ature stiffness is proportional to t\ and to (Ji — J2) 2 , 
that is, the square of the difference of electric current 
carried by a quasiparticle on each layer. We proposed 
that experiments with bilayers or heterostructures can 
serve as a new kind of probe of the effective quasiparticle 
and its doping dependence in cuprates. 
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APPENDIX A: PERTURBATION THEORY IN 
THE INTER-LAYER COUPLING 

Here we use straightforward perturbation theory of 
in the inter-layer tunneling to compute the zero temper- 
ature stiffness of the bilayer. This is an alternative to 
the effective Hamiltonian approach used in section IIVDI 
to obtain the zero temperature stiffness and provides a 
check of the results. The expansion is separated to dif- 
ferent regions in the Brillouin zone where different sets of 
levels may be nearly degenerate. For example we describe 
the expansion for wave- vectors near the Fermi surface of 
the metallic layer (layer 2). The (non-normalized) eigen- 
vectors of ([5]) corrected to first order in t±/E\ and to 
lowest order in E2/E1 are given by: 
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(Al) 



Note that the subscripts k of Ex, u, and v are suppressed 
for notational simplicity. Here u = [(1 + ^i/£'i)/2] 1 / 2 , 
v = [(1 - fc/JSi)/^, u = [(1 + i 2 /E 2 )/2} 1 /*, v = [(l- 
i 2 /E 2 )/2]y 2 , where f 2 = - ^t 2 ± /E 2 and E 2 = [§ + 
A 2 t 4 L /Ef] 1 / 2 . In addition we denote f± = 2uvv — (u 2 — 
v 2 )u and f 2 = 2uvu + (u 2 — v 2 )v. 

These states can now be used to compute the response 
to an external vector potential coupled to the Hamilto- 
nian ([5|). 

Hmf(A) = Yj\ttk + + *l>H^)*k] 

k 
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where In particular the diamagnetic response at T — in- 

^ eludes the first order correction to the ground state en- 

h k (A) = hk(0) - J A + -^KA 2 (A2) crgy in the quadratic coupling term, 

with /ik(0) given by (|6|) and 



+ & + (2|^|2) + (4|£|4 
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6k 



(A4) 
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The zero temperature paramagnetic contribution is given by the second order perturbation theory in the linear 
coupling term: 
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tj_ sin 2 k x 
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Eikie 



2k 



(A5) 



The perturbative corrections to the zero temperature 
stiffness to order r 2 = t\/E 2 k consist of two contri- 
butions. The diamagnetic contribution, proportional to 
2tr 2 (with t = ti = tz) and the paramagnetic contribu- 
tion, proportional to 2(2i) 2 r 2 / 'Ex and thus larger by a 
factor of w At/Ei ~ hp /5k (6k is the fermi surface mis- 
match of the two layers). This consideration allowed us 
to keep only the perturbative terms ~ (J\ — J2) 2 in the 
derivation of the effective Hamiltonians (jTTJ) and (fT51) . 

The zero temperature paramagnetic response scales as 
(J 1 — J2) 2 to all orders in t±. To see this we denote the 
exact set of four eigenvectors of the Hamiltonian ((6|) by 
{ I n ) } and their corresponding energies by £ n . We note 
that the matrix J can be rewritten as ( J\ — J^d + 3%! 
where / is the 4x4 unit matrix and 

*-C?S)- < A6 > 



The paramagnetic response is given by the second or- 
der correction of the ground state energy in the presence 
of an external field (as in equation (|A5|) ) . and thus in- 
volves only off-diagonal matrix elements of J in the basis 
{ \n}}. As a result, matrix elements of J2I vanish due 
to orthogonality of the eigenvectors and we are left with 



r oSTfl 7^2 V I ( n \j\ m ) I 2 / a n\ 

k neneg C " Cm 



where n € neg denotes the negative energy eigenstates. 
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